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O . ABSTRACT 
(N ; 

bJ)- Many luminous blazars which are associated with quasar-type active galac- 

^ ■ tic nuclei display broad-band spectra characterized by a large luminosity ratio 

of their high-energy (7-ray) and low-energy (synchrotron) spectral components. 
^ ■ This large ratio, reaching values up to 100, challenges the standard synchrotron 

^ self-Compton models by means of substantial departures from the minimum 

O ' power condition. Luminous blazars have also typically very hard X-ray spec- 

tra, and those in turn seem to challenge hadronic scenarios for the high energy 
• blazar emission. As shown in this paper, no such problems are faced by the mod- 

I ' els which involve Comptonization of radiation provided by a broad-line-region, 

^ ■ or dusty molecular torus. The lack or weakness of bulk Compton and Klein- 

^ ■ Nishina features indicated by the presently available data favors production of 

7-rays via up-scattering of infrared photons from hot dust. This implies that 
CN ■ the blazar emission zone is located at parsec-scale distances from the nucleus, 

and as such is possibly associated with the extended, quasi-stationary reconfine- 
ment shocks formed in relativistic outflows. This scenario predicts characteristic 
■ timescales for flux changes in luminous blazars to be days/weeks, consistent with 

T^ij- ■ the variability patterns observed in such systems at infrared, optical and 7-ray 

^ . frequencies. We also propose that the parsec-scale blazar activity can be occa- 

O ■ sionally accompanied by dissipative events taking place at sub-parsec distances 

J> . and powered by internal shocks and/or reconnection of magnetic fields. These 

^ . could account for the multiwavelength intra-day flares occasionally observed in 

. powerful blazars sources. 
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1. Introduction 

Blazars arc active galactic nuclei (AGN) dominated by Doppler boosted radiation of rel- 
ativistic jets, and as such, they provide an exceptional opportunity for studying the physics 
of innermost portions of AGN jets. However, in order to take advantage of this opportunity, 
the appropriate radiative processes responsible for generation of a broad-band jet emission 
must be identified. And, while the low-energy component of double-peaked blazar spectra is 
uniquely identified as originating via the synchrotron radiation of ultrarelativistic electrons 
(hereafter by electrons we mean both electrons and positrons) . the origin of the high-energy 
spectral component is still under debate. Possible contributions inchide inverse-Compton 
(IC) emission by electrons directly accelerated within the jet (and producing the observed 
synchrotron continuum), synchrotron radiation of pair cascades powered by hadronic pro- 
cesses, and synchrotron emission of ultra-high-energy protons and muons (see reviews by 
Sikora & Madejski 2001; Levinson 2006; Bottcher 2007). Seed photons for the IC process 
may be provided by the jet synchrotron emission as well as by the external photon sources, 
such as an accretion disk, broad line regions (BLR) and/or dusty tori (hereafter 'hot dust 
region', HDR). The IC models involving Comptonization of local synchrotron radiation are 
called synchrotron-self- Compt on (SSC) models, and those involving Comptonization of ex- 
ternal radiation - external-radiation-Compton, or, for short, external-Compton (EC) models. 

In this paper we focus on blazars with a dense radiative circumnuclear environment, 
i.e. those hosted by quasars. High-energy (HE; 7-ray) luminosities of many of those objects 
exceed their low-energy (LE; synchrotron) luminosities by a factor q = Lhe/Lle = 10 — 100. 
At the same time, their X-ray spectra are often very hard, in a number of cases character- 
ized by energy spectral indices < 0.5 (see Table 1 and references therein). Constraints 
resulting from such properties are discussed in § 2 and § 3, respectively. In addition, we 
investigate constraints imposed: by the spectral localization of the two prominent blazar 
radiative components (§4); by the lack of theoretically predicted bulk-Compton features 
in blazar spectra (§5); and by the lack of signatures of the Klein-Nishina features in their 
synchrotron continua (§6). We discuss EC(BLR) and EC(HDR) models in §7, highlighting 
their similarities and differences. We discuss the possibility of detecting electromagnetic im- 
prints of the presence of protons in the blazar jet spectra in §8. Our final results regarding 
physics of the parsec-scale jets of quasars are presented in § 9, and we briefly conclude in § 10. 
Throughout the paper the approximation 6j = Oohs = l/T is used, where dj is the opening 
angle of a conical jet, Ogbs is the angle between the line of sight and the jet axis, and V is the 
jet Lorentz factor. We also neglect the redshift dependence in the derived expressions, but 
note that this dependence goes at most as a factor (1 -|- 2;). 
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2. Lcirge Luminosity Ratio 

Since our understanding of particle acceleration processes in relativistic plasma is lim- 
ited (and relies primarily on results obtained strictly in the test-particle limit), the existing 
hadronic models are not sufficiently quantitative to make robust predictions regarding the 
relative power required to be injected in the form of high energy protons and ultrarelativistic 
electrons. Hence, the observed large luminosity ratio between the high- and the low-energy 
components in blazar spectra can be scaled by a free parameter in those models. In contrast, 
in leptonic scenarios both components are produced by the same population of ultrarelativis- 
tic electrons, and therefore the observed large values of q provide valuable constraints on the 
blazar parameters. 

Noting that for luminous, quasar-associated blazars, the Klein-Nishina and the pair 
production effects are likely to be insignificant (Moderski et al. 2005; see also §6), and that 
the effect of the anisotropy of external-Compton scattering in the source co-moving frame is 
negligible for 9,j,,^ ~ l/F (Dermer 1995; Moderski et al. 2003), the following ratios of energy 
densities can be approximated as being equal: 

'^'SSC ^ '^'syn ^ u'eC — a (-]\ 
"•syn ^ext 

where A is the 'Thomson amplification factor', u'^ is the comoving magnetic energy density, 
while u'^y^^, v'gg^, Mgxt ^-nd u'^^ are the energy densities of the synchrotron, SSC, external 
(BLR or HDR), and EC photon fields, respectively, all as measured in the source rest frame. 



2.1. Large Luminosity Ratio in the SSC Model 

In the case of the SSC emission dominating the high-energy spectral component, Lhe 
Lssc, one has 

'^SSC _ USSC _ '"'syn 



and 



Lssc u'sSC '^'syn , 

q ~ ~ ~ — J— — A, (2) 

Lsyn '^'syn ^B 

= = (3) 

'^SSC '^syn '^SSC 

where Lgyn and Lssc are the observed synchrotron and SSC luminosities, respectively. Not- 
ing next that the approximate relation between the ratio of total energy densities stored in 
non-thermal photons, w^^^, and electrons, u'^, is 

'^'rad ^ Vrad /^x 
K 1 - Vrad ' 
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where rjrad is the radiative efficiency of the injected electrons, and that in the SSC model 
with 1 one has u'^.^ ~ fJ'SSC^ '^^ obtain 

Since for the luminous blazars rjrad > 0.5, the 'high-g states' of the SSC-dominated 7- 
ray emission can be reached only for u'^/u'^ > S> 1, i.e. for conditions far away from 
equipartition. 



2.2. Lcirge Luminosity Ratio in the EC Models 

In the case of the high-energy emission dominated by the EC radiation we have Lhe — 
Lec, leading to 

^ T ~ 11' ~ 11' ' 

^syn "'syn "'B 

and ^ 

'^B ^ '^B '^syn _ _i 



"^'eC '^'syn '^'eC 



A-'q-'. (7) 



Hence, with the further approximation u'^^^ ~ u'^q (as justified for ^ 1), we obtain 

u'j^ u'j^Q 1 Tjrad 



< '^'ec < Aq {l-r]rad) 



(8) 



Obviously, in the EC models the quantities A and q are not equivalent, and the equipartition 
between magnetic field and electron energy densities can be reached for any (large) value of 
the q parameter provided 

. 1 Vrad 



? (1 - Vrad) 

Moreover, since A ~ u'^cl'^'ext^ where 



(9) 



and 



one may find that 



/ ^EC 



p2t T 

I ^ C,ext-L-'disk /i i\ 



S.ext-'-' disk'- 
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In the above, we have introduced the distance of a blazar emission zone from the nucleus 
r, the radius of the emission zone R ~ rOj ~ r/T, the accretion disk luminosity Ldisc, and 
the fraction C,ext of the disk radiation reprocessed in the BLR or HDR. This, when combined 
with Eq. ([H]), shows that the minimum power condition u'^ ~ u'^ implies bulk Lorentz factors 
of the emitting plasma 

1/4 

r~20x ( , 1 ( ,..r""'' . 1 {- — ^1 . (13) 



lO/VO.l/ VlO^^ergs-V VlO^^ergs-V V Vrad 

We note that typically in blazar sources the disk luminosity L^isk cannot be observed directly 
but instead can be estimated from the observed luminosity of broad emission lines, leading 
to Ldisk ^ 10^^ erg s~^. Also, the derived 'equipartition' value of the jet Lorentz factor, which 
depends only weakly on the model parameters, is nicely consistent with the typical values 
imphed by the observed superluminal expansion of pc/kpc-scale radio blobs in 7-ray bright, 
quasar-hosted blazars (Jorstad et al. 2001a, 2001b; Lister et al. 2009; Kovalev et al. 2009). 



3. Hard X-ray Spectra 

About 30% of blazars detected by BeppoSAX up to the energy of 50keV are character- 
ized by very hard X-ray spectra, with the X-ray spectral indices smaller than the 'canonical' 
value 0.5 (Giommi et al. 2002; Donato et al. 2005; see also the Table 1). In several cases 
strong indications were found for such hard spectra to extend up to at least ~ 100 keV. Such 
hard spectra must be produced in a slow cooling regime for the radiating particles, and below 
we investigate what constraints such observations impose on different emission models. 

3.1. Hard X-ray Spectra in the Hadronic Models 

3.1.1. Proton-Initiated Cascades 

In the hadronic models for luminous blazars, the observed X-ray emission is produced 
by synchrotron radiation of pair cascades powered by ultrarelativistic protons (Mannheim 
& Biermann 1992; Mannheim 1993). In order to obtain ax < 0.5, the cooling break in the 
electron distribution formed by such cascades, corresponding to the cooling break frequency 
in the synchrotron continuum Vsyn,c as measured in the observer frame, must correspond to 
the electron Lorentz factor 

7e ^ 6 X 10^ B'^-'" ( 1^!^] ( L] , (14) 

' ^ Viookevy V20/ ^ ' 
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where the comoving magnetic intensity is expressed in units of Gauss, B'q = B'/G. Mean- 
while, equating the time scale of the synchrotron (electron) energy losses t'^^^^ ~ 'j/\'jsyn\ 
with the dynamical time scale of the emission region t'^^,^ ~ (R/c) ~ r/(cr), one may find 

m^c^ r 

7c ^ r • (15) 

Hence, hard X-ray spectra observed in luminous blazars, if powered by hadronic interactions, 
require jet magnetic field as low as 

B'c^5x 10- (-](-) ( . (16) 

However, in such weak magnetic fields, the SSC spectral component produced by copious, 
directly accelerated electrons can lead to overproduction of the X-ray flux. In order to avoid 
this situation, the blazar emission zone must be located at sufficiently large distances to 
satisfy A = u'^y^/u'^ < q, where the comoving energy density of the the jet photons is 

- - 7 X 1°-' (lojfe) {^y i^y (1^) 

This condition, together with Eq. f|T6l) . implies in turn 

r>20( ^^^1 ( ( ^] " (^) PC . (18) 
VlO^^ergs-V VlOOkeVy V20y VlO/ ^ ^ ' 

At such distances, however, pair cascades are likely to be produced very inefficiently. 
In order to verify this, we calculate below the efficiency of the photo-meson process which 
dominates the powering of the cascades. The opacity for the photo-meson process involving 
protons with random Lorentz factor 7p is 

( \ — 1j!IL / K \ ^^P^(^P-^ (-\(X\ 

dyn 

where (ap^Kp^) ^ 0.7 x 10"^^ cm^ is the product of the photo-meson cross-section and 
inelasticity parameter averaged over the resonant energy range (Begelman, Rudak, & Sikora 
1990), n'pf^^'jp) = J^, ^ n'p^ ^idv' is the number density of photons with energies above the 
photo-meson threshold hv[i^{'-jp) ~ m-^c? j^p^ and m^c^ ~ 140 MeV is the rest-mass energy of 
a pion. For Tp^(7p) < 1 (which can be verified a posteriori) , the total efficiency of the process 
is 

/"Ypi max t 1 
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where n'^^ is the energy distribution of relativistic protons. The maximal available pro- 
ton Lorentz factor, in turn, 7p,max5 may be evaluated very roughly through the condition 
^'accilv) ~ ^'dyn (since the dynamical time scale is expected to be much shorter than any 
proton cooling time scale), where 

Cc(7p)^^^10-S/5g"'s (21) 

is the proton acceleration timescale, ~ jprnpC^ / {eB') is the proton Larmor radius, and 
the factor / accounts for the specific particle acceleration model assumed, being generally 
expected to be /> 10 (see Aharonian et al. 2002). This, together with Eq. ( |T6l) . gives 



'^p,max — OXiU IDq 



(22) 



lOy \pcj \20 J \10 J \20pcJ VlOOkeV 

and the corresponding photon threshold energy is hu[f^{;~^p^rnax) ~ (7p,max/10^)^^ eV. 



The target soft photons for the photo-meson production are provided by the synchrotron 
radiation of directly accelerated jet electrons, as well as by the external photon fields, domi- 
nated at parsec distances by hot/warm dust located in the surrounding molecular tori. The 
number density of synchrotron photons with energies larger than hvth{'lp,max) is 



n 



p, max j 



5xlO^C(7pma.) - ^s,n __ ^ 

where C(7p,max) - Lsynii" > T i>'th{.lp,max)) / Lsyn which for the typical synchrotron spectra of 
luminous blazars and 7p,m.aa; ^ 10^ is expected to be -C 1. 

As for the external target photon field, we note that the dust located at r ~ 20 pc 
and heated by the accretion disk with luminosity Lctisk ~ lO^^ergs"^ has a temperature 
T ~ 400 K (see Eq. [39] below) . Its electromagnetic spectrum peaks at huext — 0.1 eV, and its 
comoving number density is 

^ext\ ( L^isk \ f ^ \ f / hUext 



vo.iy vio^^ergs-v v2oy v20pcy vo-ievy ^ ^ 

(see Eq. HOl below). Hence, the target radiation field for the photo-meson process is dominated 
by the external radiation, and thus the opacity for protons with energies 7p > 'jp^th — 
itlt^c^ / {T hvext) can be written as 



(f ) - 2 X 10- ( { {^Y{ " ■ (25) 



0.1 y VlO^'^ergs-V V20pc/ 
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For the most optimistic acceleration scenario with / = 10, and for the proton energy distri- 
bution of the form n'^^ oc 7"^, this gives the total efficiency 



7?^f ) ~ tI,T^ ^^fa''"""/^^'^^^ ^ 2 X 10-^ . (26) 
111 7p, 



With such a low efficiency of the photo-meson production process, the jet power required to 
explain the observed 7-ray luminosities of the order of 10*^^ergs~^ is 

where rjdiss is the fraction of the jet energy flux dissipated in the blazar emission zone, and 
r)p is the fraction of the dissipated energy channeled into relativistic protons. The required 
jet kinetic luminosity is therefore orders of magnitude larger than the Eddington luminos- 
ity of a black hole with mass 10^ Mq. Such serious energetic difficulties could be overcome 
only by postulating ad hoc some substructure of the blazar emission zone, involving spa- 
tial separation, and very different physical conditions characterizing regions where hadronic 
and leptonic processes operate. However, several arbitrary assumptions involved in such 
a scenario as well as an increased number of free parameters would make this much less 
attractive. 



3.1.2. Proton- Synchrotron Emission 

Alternatively, in order to explain hard X-ray spectra of luminous blazars still in the 
framework of hadronic models, one may postulate that the entire observed high-energy (X- 
ray-to-7-ray) flux is dominated by the direct synchrotron radiation of ultrarelativistic pro- 
tons. Noting the relation between the characteristic proton- and electron-related synchrotron 
frequencies Vp^synilp = 7e) = (iTie/mp) Ue^syn, and the relation between the corresponding 
cooling rates \'yp\syn{lp = 7e) = (^e/'^^p)^ |7e|sj(n (see, e.g., Aharonian 2000), one can find 
the appropriate radiative efficiency around the 7-ray luminosity peak 

...,„(7„„.) - % - 2 . (^) . (28) 



with 



^^,^,:,3 y.UfB'^-"'(^y' . (29) 



The resulting efficiency of the synchrotron proton emission is, again, very low, and can be 
increased only by assuming very large magnetic field intensity within the emission zone. 
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However, since the magnetic energy flux is limited by the total jet kinetic power, Lb — 
cu'^Tir'^ < Lj, magnetic fields Bq S> 1 are expected only at r <^ Ipc. In this context, one 
has to keep in mind that at small distances from the black hole, AGN jets are still in the 
acceleration phase, and thus the radiation produced there is not boosted sufficiently strongly 
to account for the observed large luminosities. The distance scale of the jet acceleration, 
which corresponds roughly with the distance of the conversion of the Poynting ffux dominated 
jet to that dominated by matter, is expected to be Vq > 10^ Rg ~ 0.03 (M^h/IO^ Mq) pc 
(Komissarov et al. 2007), where Rg is the radius of the central black hole. Considering those 
additional constraints, we obtain 

/ . -1/2 X -p X -3/2 /, \V2 / T \3/4 

X,„„(W) < 2 X 10-3 (^) (r^j ^_L^^ _j3^^ 

where Lj ~ 10^'' erg s^^ corresponds to the Eddington luminosity of the 10^ Mq black hole. 
With such an efficiency, the proton synchrotron peak luminosity is 

Lp ,,„(10MeV) < 4 X 10^^ ( ^) ( -^^i^] ( iV ( —h \ ergs'^ , 

^''^"^ ' V 0.1 / Vo.5/ V2X10-V V20y VlO'^ ergs-V 

(31) 

i.e. orders of magnitude below the observed 7-ray luminosity of powerful blazars. 



3.2. Hard X-ray Spectra in the SSC Model 

In the SSC model for large-g blazars, the electron coohng rate is dominated by Comp- 
tonization of co-spatially produced synchrotron radiation, and therefore reads as 

l7l-l7|5SC- , (32) 

where 

Balancing next the corresponding comoving time scale of radiative losses t'^^^^ ~ 'y/\'j\ssc 
with the dynamical time scale of the emission region t'^y^, one may find the critical electron 
energy corresponding to the transition between fast- and slow-cooling regimes. 



Electrons with 7 < 7c cool inefficiently, Compton-upscattering the far-infrared synchrotron 
radiation up to MeV photon energies. Therefore, the observed hard X-ray spectra can be 
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reproduced by the SSC model in a framework of the standard approach involving injection of 
some given (single or broken) power-law electron energy distribution to the emission region 
and its subsequent radiative cooling, provided the low-energy segment of the injected electron 
energy distribution is not softer that oc 7"^. 



3.3. Hard X-ray Spectra in the EC Models 

In the EC models for blazars with large values of the q parameter, radiative cooling 
of electrons is dominated by Comptonization of the external radiation fields. As before, by 
equating the appropriate cooling time scale 

t'cooi ^ ^J, 35 

with the dynamical time scale of the emission region t^^^^, one finds the critical electron 
Lorentz factor 

Ic -, — , 36 

where we assumed that the energy density of the external photon field may be approximated 
as being constant up to some characteristic distance Text-, and that the blazar emission zone 
is located at ro < r < Vexti where tq is the distance scale of the jet acceleration (see § 3.1.2). 

In the case of the dominant EC(BLR) emission, i.e. when the electron cooling is dom- 
inated by Comptonization of the photon field due to the broad emission lines, we have 
ro < r < Text = tblr, where 

(Ghisellini & Tavecchio 2008 and references therein), and 



/ -p2 iBLRLdisk^^ TO f ^BLr\ f ^ \^ -3 /oo\ 



BLR 

If, instead, tblr < r < thdr, where 



^""^"'(. lO^^eTgs-O iT^J ^''^ 



the electron cooling is dominated by IC up-scattering of the near-IR photons emitted by 
the hot dust with the temperature T, located in the molecular torus and irradiated by the 
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accretion disk (Nenkova et al. 2008 and references therein). In such a case 



u 



UHDR^ ~ ——2 — ~ 8 X 10 M — — — — — erg cm ^ (40) 



e^t--'MUH^ - 4T,rljyj^c VIO^K; V 0.1 ; V20/ 

and the EC(HDR) emission dominates the high-energy radiative output of the source. 
With all the values quoted above, the critical electron Lorentz factor reads as 



20 J VlO^'ergs-V )^ 250 (r/r^Di?)"' (CWO-l)"' (T/IO^K)"'-" 

and the corresponding break frequency in the high-energy spectral component, j^ec, 
'j^T^Uext, may be evaluated as 

hvEc,c\ ( Ldisk J 0.2 (r/rBLfi)"^ (^blr/0.1' 



(41) 



MeV ; - VlO^^ergs-V \ 8 {r/rHDRy' {^hdr/O-I)-' {T/lO'Ky'-' ' ^ ' 

where hugxt has been substituted with hvBLR — 10 eV and hvHDR — 0.3 eV, respectively. 
Hence, the production of high-energy emission in the EC(HDR) model proceeds in the slow 
cooling regime up to the MeV photon energy range, in agreement with the position of 
luminosity peaks observed in luminous blazars. In the EC(BLR) model the 7-ray spectra 
are expected to soften below the MeV range, but still at sufficiently high photon energies to 
assure consistency with observations. Hence, in similarity to the case of the SSC model, the 
observed very hard X-ray spectra of powerful blazar sources are in agreement with the EC 
scenario for the high-energy jet emission. 

High energy radiation produced in the innermost portions of AGN jets may also be 
contributed by Comptonization of the direct emission of the accretion disk illuminating the 
blazar emission zone from the sides and from behind. This contribution should in particular 
dominate the high-energy radiative output of the jet at distances 

. < r.„.,«.. 0.6 (^^— j =.0.01 (^—j (yoS^J \wm;J 

(43) 

(Dermer & Schlickeiser 2002). Energy dissipation events taking place on such small scales 
may be responsible for production of intraday flares observed sporadically in some blazar 
sources (see § 7.2). 



4. Synchrotron and Inverse-Compton Peak Frequencies 



In the case of luminous blazars, the observed synchrotron peak frequencies seem to 
cluster within a relatively narrow frequency range 10^^ — 10^"^ Hz (Ghisellini & Tavecchio 
2008 and references therein). Meanwhile, the peak of the high-energy spectral components of 
the objects discussed here have to be located between 0.1 MeV and 0.1 GeV photon energies, 
as indicated by the observed values of the X-ray and 7-ray spectral indices (e.g., Fossati et 
al. 1998, Ghisellini & Tavecchio 2008). Again, the relative location of the low- and high- 
energy spectral peaks does not provide any constraints on the hadronic models, because the 
poorly known particle acceleration efficiency (the / parameter introduced in Eq. HT] above) is 
claimed to be substantially different for electrons and protons (e.g., Reimer et al. 2004), and 
therefore the ratio of peak frequencies is in fact a free parameter. However, in the leptonic 
scenarios for the high-energy blazar emission, the relative location of spectral peaks offers 
interesting constraints on the model parameters, as discussed below. 



4.1. Peak Frequencies in the SSC Model 

Peak frequencies of the synchrotron and SSC components in the spectra of blazar sources 
are determined by the critical (break) energies of electrons emitting bulk of the emission. 



7..-,/^^. (44) 

y ^syn, br 

When combined with the standard formulae for the corresponding break frequency of the 
synchrotron photons 

^4 X 10*^5^72^.1 Hz, (45) 

we obtain 

In the SSC model, magnetic field intensity B' given by the relation u'^y^/u'^ ~ g is 



10^7 erg s-V V20/ \pc J VlO 
and hence 

This would then be consistent with the observational constraints for the location of the 
emitting zone r > a few pc. 
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However, the SSC mechanism can dominate production of the 7-ray emission only if 
the energy density of synchrotron radiation is larger than energy density of the external 
radiation field, both as measured in the source co-moving frame, i.e. if u'^y^ > u'^^^ ~ u^^t^'^- 
This condition can be satisfied on parsec scales, provided the contribution from the BLR 
at these distances is already negligible (as expected even in luminous blazars) and the dust 
temperature is much below the sublimation temperature. The second requirement can be 
verified observationally by investigating the presence of a hot dust in mid- and near-IR data 
for non-blazar radio-loud quasars, i.e. for jet sources observed at 9obs ^ 1/r. According to 
Cleary et al. (2007), who analyzed Spitzer observations of 3C quasars, the dust temperatures 
are in the range T ~ 100 — 300 K. Such temperatures are derived assuming that the observed 
mid-IR excess in the spectra of analyzed objects is produced by a high energy tail of the 
synchrotron radiation of a large scale jet. If this is indeed the case, then there is a distance 
range between Tblr and Thdr where the SSC process can dominate the production of 7-rays 
over the EC(HDR) mechanism. On the other hand, we note that the near-lR quasar continua 
join the harder optical ones almost always around the same frequency, namely ~ 1 /xm, in 
both radio-quiet and radio-loud sources (Barvainis 1987; Ghkman et al. 2006; Netzer et al. 
2007; Labita, Treves, & Falomo 2008). Therefore, the assumption regarding the synchrotron 
origin of near IR radiation seems not to be justified, requiring in addition very precise tuning 
of the jet parameters. At the same time, a natural explanation for the observed spectral 
dip around 1 frequencies in quasar spectra is provided by the clumpy torus model, in a 
framework of which the near-lR emission is produced by a hot dust with the temperature 
close to the sublimation value, T ~ 1500 K (Nenkova et al. 2008). 



4.2. Peeik; Frequencies in the EC Models 

In the EC models, the high energy (inverse-Compton) spectral component is expected 
to peak at the observed frequencies 

J^ECbr - Ibr^'^^ext , (49) 

and the ratio of the peak frequencies reads as 

^ 10^ ( -] ( ^] B'^-' . (50) 

With the magnetic field given by the q — T^Uext/u'^ condition, 

8^ ( u,,t / r \ / g X -1/2 f 5 (Cblr/O-I)'/' 



B'r 




^ '\l q Vergcm-3; ^\2Q)\iq) \ Q.l {iHDR/^.lf'^T /WKf 



(51) 
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this ratio is then surprisingly similar for both cases considered here, EC(BLR) and EC(HDR): 

i^synM ~ ^10^ \ (CWO.l)-^/^ (T/IO^K)-^-^ ' ^^^^ 

leading to the expected spectral location of the EC luminosity peak 

/...a.-0.l(JL)"^(|2)"'^(ii^)GeV, (63) 

with the assumed T ~ 10'^ K. The hvEC,br evaluated above seems to be somewhat higher than 
the peak photon energies of the high-energy continua claimed for luminous blazars (GhiseUini 
& Tavecchio 2008), as long as VsynM- is not smaller than 10^^ Hz. However, we note that 
the high energy peak frequencies are not observed directly, but only reconstructed from the 
power-law slopes of the lower-energy (X-ray) and higher-energy (GeV) segments of the IC 
continua, assuming a single-broken power-law spectral shape of the emission. Meanwhile, if 
the injected electron spectrum is flatter than below electron energies 75^ corresponding 
to the peak frequencies Vsyn^br and VECbr (as seems to be indicated by the observed hard 
X-ray spectra with ax < 0.5), and steeper than that above 7^^ (see in this context GhiseUini 
& Tavecchio 2008, but also §9 of this paper), the 'reconstructed' peak frequency can easily 
be lower than the 'true' break frequency vecm- ^ factor of (^'Ec,c/^EC,6r)'-°'^~""'^''^"^~"'^\ 
i.e. by possibly as much as ~ 10 for the typical values of the X-ray and GeV spectral indices 
of luminous blazars (see Table 1). 



5. Lack of Bulk-Compton Spectral Features 

In astrophysical jets, all ultrarelativistic particles cool due to radiative energy losses and 
also due to adiabatic expansion of the emitting regions. Hence, the observed high-energy 
non-thermal radiation requires in situ particle acceleration and one may expect to observe 
at any given instant radiation from both 'cold' and 'hot' population of particles. Typically, 
jet electrons constituting a cold population are considered to be non-relativistic in the jet 
rest frame in the innermost parts of the outflow (say, hundreds of Schwarzschild radii from 
the active center), because of the extremely rapid (catastrophic) particle cooling at the jet 
base. Such electrons, if indeed present in a steady flow, are predicted to produce distinct 
bulk-Compton features in blazar spectra, most pronounced around hui^c ~ V'^hv^xt photon 
energies, i.e. around ~ 1 — 4keV for F ~ 10 — 20 and hvext ~ lOeV (Begelman & Sikora 
1987; Sikora et al. 1997; Sikora & Madejski 2000). In addition, in the internal shock model 
(Spada et al. 2001), presence of cold jet electrons should manifest itself as soft and hard X- 
ray precursors of non-thermal blazar flares (Moderski et al. 2004). Both those 'steady' and 
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'dynamic' bulk-Compton features are expected to be particularly strong if produced close to 
the black hole, in the dense radiation field of an accretion disk. So far, no clear detection 
of such features has been reported, with the possible exception of weak soft X-ray excesses 
claimed in the spectra of some blazars (see, e.g., Worsley et al. 2004; Kataoka et al. 2008; De 
Rosa et al. 2008). Such an apparent weakness of bulk-Compton features may result from the 
fact that in the very central regions of AGN jet the (magnetohydrodynamical) acceleration 
of the outflow is not yet fully complete (as noted previously in Sikora et al. 2005; Celotti, 
Ghisellini & Fabian 2007). On the other hand, analogous spectral features produced further 
out at the point where the jet is already accelerated to terminal bulk velocities are predicted 
to be still detectable, in particular if the blazar emission zone is enclosed within the BLR. 

However, it should be emphasized that the involved assumptions regarding cold electron 
population may not be valid at all scales. Namely, at some larger distances from the nu- 
cleus, the radiative and adiabatic cooling are expected to become inefficient when compared 
to the acceleration processes at very low electron energies. As a result, shock heating and/or 
interaction of particles with magnetic turbulence (driven inevitably by MHD instabilities 
of a magnetized outflow, and the related magnetic reconnection processes) may effectively 
accelerate all the cold electrons to energies 7^1, replacing the non-relativistic Maxwellian- 
like electron population with a very hard, quasi-power-law distribution that joins smoothly 
the higher-energy electron tail (see §9 below). In such a case, bulk-Compton features are 
expected to be absent. Yet another alternative explanation for the apparent lack of radia- 
tive signatures of cold electrons in blazar spectra is that formation and acceleration of jets 
in luminous sources is still not completed within the BLR. If so, bulk-Compton features 
produced at r < tblr may be simply too weak to be distinguished from the broad-band 
continuum generated via EC(HDR) process. The latter possibility is supported by the fact 
that in several cases pure power-law X-ray continua extend down to energies 0.1 — 0.3 keV, 
with no evidence for any soft excess (e.g. in 3C 279: see Lawson & M'^Hardy 1998; CoUmar 
et al. 2007). 



6. Lack of Klein-Nishina Spectral Features 

Recent studies of the Klein-Nishina (KN) effects of relativistic sources immersed in a 
dense external radiation fleld and dominated by the IC emission show that for the maximum 
electron energies 'jmax 3> 'Jkn = n^eC^ /4:ThUc.j.f the IC continuum should break at much larger 
frequencies than ukn = IxN^^^ext, and that the synchrotron component should harden 
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around 

u,y^,KN^^xlO'l'KM^B'QRzc^3xlO''B'Q (^^^ (^^^ Hz. (54) 

(Dcrmcr & Atoyan 2002; Moderski et al. 2005). Both these effects result from hardening of 
the injected electron energy distribution due to the reduction of the IC cooling rate in the 
KN regime, i.e. for the electrons with 7 > '^kn- In a framework of the EC(BLR) model, 
the related spectral features are expected to be pronounced already for electrons with ener- 
gies 7 > 3 X 10^, and the hardening of the synchrotron continuum (being especially strong 
for mono-energetic external radiation fields) is predicted to be visible already in the optical 
band. However, as shown by Tavecchio & Ghisellini (2008), emission of BLR is not exactly 
monochromatic, and in particular is characterized by a significant low-frequency power-law 
tail. As a result, the expected KN hardening in the optical band may be relatively weak, 
though possibly still prominent in the UV band. The dependence of such KN effects on 
the maximum electron energy and the parameter q for the EC (BLR) model is illustrated in 
Figures 1 and 2, with an assumed electron energy distribution A^^ oc 7"* with s — 2.4. As 
shown in these Figures, high-g blazars with synchrotron and IC spectra extending respec- 
tively above 1 eV and 10 GeV photon energies should display in their broad-band continua 
the KN spectral features discussed above, unless the production of 7-rays is dominated by 
Comptonization of IR radiation of hot dust, and only if the electron energy distribution 
extends to 7 > 3 x 10^ with the power-law slope not much steeper than s ~ 2.5. This 
provides a possibility to determine the characteristic energies of the dominant seed photons 
via analysis of the Klein-Nishina effects, manifested by different spectral shapes of the high 
energy tails of the synchrotron and EC continua. 



7. EC(HDR) versus EC (BLR) Models 

As shown in §4, the same luminosity ratio of the IC and synchrotron spectral compo- 
nents implies similar high-energy peak frequencies in both EC(BLR) and EC(HDR) models. 
And this is clearly coincidental since the location of the two spectral peaks happens to scale 
with a factor Vextl (see §3.3). Interestingly, the conspiracy between the two models is 
reflected also by several other - though not by all - important model output parameters, as 
discussed below. 
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7.1. Similcirities 

Let us consider a 'fiducial' source located within a given distance range Ar ~ r, and 
immersed in a given external radiation field. The observed X-ray and 7-ray spectral indices 
of such a source, and respectively, allow one to reconstruct the energy distribution 
of radiating electrons, approximated here for simplicity by a single-broken power-law. In 
particular, using the approximate formula 



L^diy ~ N^d'y I7I m^c^ (55) 



we find that 



7.. = {^] , (57) 



N, = CN,l-''{l+lhbry'-'' , (56) 

'^x '^obs,x 

Si = 2ax + 1 , and S2 = 2a^ + 1 . (62) 

Here F,y^^^ are the observed monochromatic energy fluxes at some particular observed fre- 
quency Uobs, dL is the luminosity distance of the source, and z is its redshift. 

The number of electrons contributing to the radiation at a given instant of observation 
can be approximated by 

K^U^d^^^, (63) 
as long as 7;,^ 3> 1 and 7^4^ = 1. Hence, the electron number flux is 

[R/TC] [rext/c) Uext Text 

and the ratio of electron energy fluxes in the two considered models reads as 



Ne, HDR 
Ne, BLR 



( YhDrV ( UhDr\ ^ ( ThdiA ( T \ / iHDR\ ^ 
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which, for the sublimation temperature T ~ 1500K and $,hdr ~ '^^blr-, is of the order of 
unity. 

Assuming further that the energy flux of a jet is dominated by nonrelativistic protons, 
the proton number flux is simply 

and therefore the pair content in the two models is 

{ne/np)HDR _ Ne,HDR 



BLR N^BLR 



(67) 



i.e. almost exactly the same. In addition, similar value for iVg derived for the two models 
implies also similar electron energy fluxes, 

Le = (7) mec' VN, , (68) 

where (7) = J N^'-fd'-f/ J N^d'-f ~ i'n.'ybr is the average electron Lorentz factor (assuming 
— 0.5), and 'jbr = ^^^^/i^ec, br/ ^ext- It follows directly from the above discussion (§3 and 
§4) that (7), and thus Lg, is the same for the EC(BLR) and EC(HDR) models. Finally, 
also magnetic energy fluxes Lb — cu'^F'^ttB? ~ cu^vrr^ do not differ in both cases considered 
above, namely 

LbJIDR ^ UHDR rluR ^ Q g3 / T\ iHDR 

Lb, BLR UblrTblr \10-^KJ ^blr 



7.2. Differences 

Certainly there are some aspects which differentiate the two models considered. These 
are (i) the critical synchrotron self-absorption frequencies, (ii) the characteristic variability 
time scales, and (iii) the low energy breaks/cut-offs of the IC spectral component. 

(i) The synchrotron self- absorption makes the source optically thick below the frequency 

u^c^Ax 10- 5^^/^ ( ( -) ( -) Hz , (70) 
^ VlO^^ergs-V V20y \pc J ' ^ ^ 



which may be re-written as 



W^m) VlO^^ergs-V \ Q.2{B'/{).lGf" {r/rHDRV^'^ 
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Meanwhile, millimeter and submillimeter observations indicate absorption break in blazar 
spectra around 1mm (see, e.g., Gear et al. 1994). This favors the EC(HDR) model over the 
EC (BLR) version. 

(ii) The shortest time scales of flares differ in the two EC models by a factor thdr/tblr ~ 
30(T/10^K)2 ^ In particular, one has 

(trr,,n\ _ ( Lg^^k \ f ^ j 0.4 



\dciy J ' \lO*^ergs-^ J \20j ^ | 10 (T/IO^ K)"''^ ' ^^^^ 

Hence, the characteristic variability time scale predicted by the EC(HDR) model is more 
consistent with the one implied by the observed variability patters in 'optically violent vari- 
able' quasars, which is of the order of days/weeks (Pica et al. 1988) A nice example where 
such time scales are seen both in the 7-ray and IR/optical frequency ranges is provided by 
multiwavelength observations of blazar 3C 454.3 (see Fig. 1 in Bonning et al. 2009). 

We note in this context that in some cases rapid, day/intraday flares are imposed on 
the observed light-curves of blazar sources. One example is the extremely rapid (~ 2 hours) 
optical flare, as recorded also in 3C 454.3 (Raiteri et al. 2008), or intraday 7-ray flare 
recorded by Fermi/L AT in the radio quasar PKS 1454-354 (Abdo et al. 2009a). Such rapid 
flares may indicate that in addition to the main 'energy dissipation site' located at parsec- 
scale distances from the active center, there are also episodic dissipation events taking place 
on sub-parsec distances, possibly related to MHD instabilities and/or internal shocks within 
the inner parts of a strongly magnetized outflow. 

(iii) Because the spectra produced by relativistic electrons emerge at energies huiow — 
T'^i'ext: in the case of the EC(BLR) model a low energy break/cut-off should be observed 
around photon energies ~ 4 (r/20)^keV, while in the case of the EC(HDR) should be at 
around ~ 0.1 (r/20)^ eV (i.e., basically below the low-energy threshold of the available X-ray 
instruments). Noting however that the X-rays can be contributed by several other processes 
(such as SSC emission, bulk-Compton radiation, or high-energy emission of accretion disk 
corona), the apparent lack of the expected break can be used as an argument in favor of the 
EC(HDR) model only with caution. 



8. Hadrons in the EC Models 

If the EC models are correct, then the broad-band blazar spectra allow us to estimate the 
number of jet electrons down to lowest energies, and to flnd the corresponding particle flux. 
However, since (by deflnition) in the 'leptonic models' the entire observed radiative output 
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of a jet is produced by primary (i.e., directly accelerated) electrons, the proton content of 
an outflow remains in principle unconstrained. The most recent analysis indicate that it 
cannot be negligible, though. In fact, several different authors came to the conclusion that 
the cold protons indeed dominate the jet bulk kinetic power, at least at > tblr scales, 
and that the number ratio of electron-proton to electron-positron pairs in quasar jets is of 
the order of 0.1 - 1 (Sikora & Madejski 2000, Celotti & Ghisellini 2008). Uncertainties 
regarding the proton content could be solved directly provided acceleration of protons up to 
ultra-relativistic energies indeed takes place in the innermost parts of quasar jets with high 
efficiency, and is followed by production of some distinctive electromagnetic features, such 
as hardening of 7-ray spectra above 10 GeV or softening of X-ray spectra at lower energies. 
Again, lack of such spectral features does not exclude the dynamical role of protons in 
general, but only limits the number of ultrarelativistic hadrons, and thus constrains the 
particle acceleration models rather than the jet content. 

In order to discuss this issue in more detail, we estimate the efficiency of the photo- 
meson production process for distances and magnetic fields appropriate for the EC models 
of luminous blazars. With the target photons provided by the local synchrotron emission of 
directly accelerated electrons, we have the appropriate opacity 



n-isyr^)(^ ) ^ (a- K ) 



Lsynlp dip) 



(73) 



where Cilp) < 1 depends on the shape of the synchrotron spectrum (see §3.1). Assuming 
further that protons are injected with the same power-law index as directly accelerated 
jet electrons producing the observed 7-ray emission in the fast cooling regime, we have 
Sp = 2a-f > 2 for the typical 7-ray spectra of luminous blazars > 1 (Abdo et al. 2009b). 
Hence, because in addition ri^'^^\sp>2 > ?7p7^"'^|sp=2 and C(7p) < 1, we find using Eq. fl20l) 
that 
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The particular value of the magnetic energy flux Lb — 10^^ ergs ^ anticipated above, is 
justified by Eqs. fl38|l . fHOl) . and fl52l) . which in turn give 
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As shown above, in the EC models for luminous blazars, high efficiency of the photo-meson 
production is expected only at distances r ^ Ipc. This implies that the maximal hadronic 
radiative output is reachable at around r ~ tq, where the jet is already fully developed (i.e., 
accelerated and collimated), and the Doppler boosting is large. This efficiency can approach 
~ (//10)~^(ro/0.03 pc)~^%. Hence, even in this most optimistic case, the efficiency of the 
photo-meson production is much lower than the ~ 50% radiative efficiency of jet electrons. 

Let us further evaluate the number densities of the external (BLR and HDR) photon 
fields in the jet rest frame, approximating these by a monoenergetic distributions with the 
characteristic comoving photon energies hv'^j^^ ~ 200(r/20)eV and hv'^jj^j^ ~ 6(r/20)eV, 
respectively: 
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Hence, this implies the total photo-meson production efficiency to be 
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for ro < r < rsLij, and 
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Summarizing this section, we conclude that some imprints of hadronic activity in the 
EC-dominated high-energy blazar spectra are anticipated only if the blazar zone is located at 
ro < r < tblr- Still, this, would require almost unrealistically efficient proton acceleration, 
with / < 10. Any hadronic activity is expected to be completely absent in the blazar 
spectra if the blazar emission zone is located at distances r > tblr (i.e. at distances where 
the EC(HDR) process dominates), or if / > 10. 



9. Discussion 

According to the analysis presented in the previous sections, the broad-band emission 
of luminous blazars cannot be easily explained in the framework of the hadronic scenarios 
or the SSC model, but instead is consistent with the EC models. In addition, the EC(HDR) 
model is favored over the EC(BLR) variant. This, in turn, fixes the location of the blazar 
emission zone at r ~ 1 — lOpc from the nucleus. Note that in the EC(BLR) model, the high 
energy blazar emission is produced deeply within the millimeter photosphere and therefore 
one should not expect simultaneous high amplitude 7-ray and mm-band variability. In 
contrast, in the EC(HDR) model the blazar zone extends up to the millimeter photosphere, 
and thus mm-band variations accompanying 7-ray flares are expected, as in fact observed, 
e.g., during the 2006 outburst event in 3C 454.3 (Krichbaum et al. 2008). Hence, systematic 
monitoring of luminous blazars at mm/submm wavelengths and GeV photon energies is of 
primary interest, allowing to confirm our conclusions. 

Regardless of the observational verification, advocating the EC(HDR) model requires 
addressing the following questions: (i) why is the blazar emission zone located at some 
particular (and, quite large: r ~ IQ^Rg for lO^M© black hole!) distance from the nucleus? 
(ii) what is the dominant particle acceleration process involved, and why is it restricted to the 
relatively narrow distance range along the outflow? These questions are in fact related, since 
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the location of the blazar emission zone is primarily determined by the energy dissipation 
processes rather than by some other factors, such as a presence of seed photons for IC 
scattering, i.e. the enhanced radiative cooling of jet electrons. To justify this statement, 
in Figure 3 we have plotted energy densities of different photon fields (and also of the jet 
magnetic field), all as measured in the jet rest frame, as functions of the distance from the 
center r. Here we assumed a fully developed outfiow with the main parameters corresponding 
to the analysis presented in the previous sections. As shown, the total energy density of 
different fields contributing to radiative cooling of jet particles, u[g^ ~ '^b^'^blr'^'^hdr'^ 
decreases along the jet roughly as oc r~^, being initially dominated by the magnetic field 
energy density (r < 0.03 pc), later by the energy density of the BLR (r ~ 0.03 — 1 pc), and 
next by the energy density of the HDR (r > Ipc). Thus, no particular distance scale is 
favored as long as solely the 'u!^g^ parameter is considered. 

The location of the blazar emission zone must be then determined by the energy dissi- 
pation/particle acceleration processes within relativistic outfiow, responsible for generation 
of the non-thermal electron population. In general (though most likely also quite naively), 
one can expect that the dominant acceleration mechanism should be related to the mag- 
netic reconnection process in the case of a jet dominated by the magnetic field, and to the 
Fermi-type processes in the case of a matter-dominated outfiow. In this context, we note 
that the most recent models for the jet formation indicate that AGN jets are in fact dom- 
inated by the Poynting fiux (and are stable against current-driven Z-pinch instabilities) up 
to at least r ~ V^^Rg (see, e.g., McKinney & Blandford 2009 and references therein). At 
further distances from the jet base, AGN outflows may be considered as fully developed 
(i.e., accelerated and coUimated to the terminal values of F and Qj parameters), and likely 
(though not necessarily) converted to matter-dominated structures due to the development 
of various (mainly kink and pinch) current-driven instabilities. Such instabilities provide a 
sink for the jet magnetic field, and energize jet particles by forming shocks within, and/or 
injecting turbulence into the outfiow. 

Not much is known about the magnetic reconnection in the relativistic regime; in partic- 
ular the accompanying particle acceleration processes remain elusive (see Lyubarsky 2005; 
Lyubarsky & Livcrts 2008; Lyutikov & Uzdensky 2003). The results of the analysis pre- 
sented in this paper, which suggest that bulk of the blazar emission is produced far away 
from the jet base, may thus be taken as an indication for a low efficiency of the electron ac- 
celeration within the innermost, magnetically-dominated parts of quasar jets. That is to say, 
the reconnection-driven acceleration processes may still operate at distances r < 10^i?g, but 
the resulting radiative output does not seem to be sufficiently strong to dominate over the 
emission produced further out from the nuclei of such luminous blazars. In other words, the 
efficiency of particle acceleration has to increase significantly at some larger distance from 
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the core. And indeed, as mentioned above, initially suppressed current-driven instabilities 
are supposed to develop starting from r ~ lO^Rg, i.e., r ~ 0.03 pc for a lO^M© black hole, 
where the BLR dominates the radiative cooling of the jet electrons. Still, this is about two 
orders of magnitude below the location of a blazar emission zone favored by our analysis. 

One can propose some ad hoc explanation for the apparent discrepancy indicated above, 
such as, for example, the stabilizing role of a velocity shear within relativistic outflow (cf. 
Mizuno et al. 2007), precluding formation of shocks or strong turbulence up to the distances 
of ~ lO'^Rg and beyond. The other possible explanation would be to postulate that the 
energy dissipation processes associated with the internal shocks formed due to steepening 
of kink instabilities around 10'^ Rg, are much less efficient in accelerating jet electrons to 
ultrarelativistic energies than the analogous mechanism associated with, e.g., extended re- 
confinement shocks formed further out from the nucleus. In fact, a role of reconfinement 
shocks in shaping radiative and morphological properties of blazar jets has been discussed 
previously by several authors (Jorstad 2001a, b; Cheung et al. 2007; Sikora et al. 2008; 
Nalewajko & Sikora 2009; Bromberg & Levinson 2009). The reason for such a difference 
in the acceleration efficiency may be due to the different turbulence/magnetic field condi- 
tions for both kink-driven and reconfinement shocks: note that the most recent analysis of 
the shock acceleration in a relativistic regime clearly shows the importance of background 
(plasma) conditions around the shock front in shaping the spectral properties of non-thermal 
particles generated there (e.g., Niemiec & Ostrowski 2004; Sironi & Spitkovsky 2009). 

But should the dominant particle acceleration process in quasar jets be undoubtedly 
identified with shocks? In fact, there are strong indications for this. Note, first, that the 
maximum energy of ultrarelativistic electrons responsible for the electromagnetic emission 
of luminous blazars is not high (7 < 10^). Meanwhile, by equating the radiative cooling 
timescale t'^^^i ~ niec/aT^ u'^^^ with the electron acceleration timescale t'^^^ ~ f ^niec/eB' , 
one gets the maximum available electron energy 7 ~ 10^ for the acceleration efficiency 
parameter as low as ~ 10~^ — 10~* at any distance 0.01 pc< r < 10 pc. In other words, 
the maximum electron energy, if limited solely by the radiative energy losses, might be 
expected to be much higher than observed, since the often invoked 'maximum' efficiency of 
particle acceleration in cosmic plasma corresponds to /'^ ~ 0.1 (see also in this context Inoue 
& Takahara 1996). Second, the electron energy distribution reconstructed from the observed 
broad-band electromagnetic spectra of luminous blazars seems to be best approximated by 
a broken power-law oc 7"'' with s < 2 for 7 < 7^ < 10^, and s > 2 for 7 > 7^. 
This manifests itself in the so-called 'blazar sequence' (Fossati et al. 1998), which, even 
though it is still pending the observational confirmation, reflects the robust finding that 
the peak electron energy in powerful, quasar-hosted blazars is always jcr ^ 10^ (see Celotti 
& Ghisellini 2008). In fact, exactly this kind of electron spectrum is expected to form at 
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relativistic (perpendicular) shock mediated by cold protons, for which the proton inertia (and 
not the radiative cooling rate!) determines the critical electron energy ^ mp/rrif. ~ 10^, 
as discussed in Stawarz et al. (2007) for the case of terminal shocks in quasar jets. 

In other words, the electron energy distribution 'reconstructed' from the broad-band 
emission spectra of luminous blazars indicates that the appropriate acceleration processes 
are related to relativistic shocks formed in a matter (proton) dominated outflow. This assures 
self-consistency of the external Compton models favored by our analysis. We note in this 
context, that the low-energy segment of the electron energy distribution discussed above (7 < 
7cr) is expected to form not due to the diffusive (Ist-ordcr Fermi) shock acceleration (which 
may operate only at 7 > 7cr), but due to interactions of relativistic electrons (with gyroradii 
smaller than Larmor radii of cold protons) with turbulence and electromagnetic waves within 
the shock front. The best studied process of this type is a resonant absorption of proton 
cyclotron emission (produced by cold hadrons reflected from a magnetized perpendicular 
shock front) by the jet electrons, as discussed first by Hoshino et al. (1992). Interestingly, 
the most recent analysis of this mechanism (Amato & Arons 2006) reveals a power-law form 
of the electron energy distribution iV^ cx 7"'' with 1 < s < 2 (depending on the plasma 
composition) within the electron energy range 1 < 7 < F^/j mp/mg, where F^/j is the shock 
bulk Lorentz factor in the upstream plasma rest frame. This is in fact in a very good 
agreement with the hard X-ray spectra observed in luminous blazars, which in turn often 
exhibit X-ray spectral indices < 0.5 (see Table 1), and with the general form of the 
electron spectrum invoked in modehng blazar spectra (Celotti & Ghisehni 2008; GhiseUini 
& Tavecchio 2008). 



10. Conclusions 

• Hadronic models cannot reproduce the very hard X-ray spectra observed in a number of 
luminous blazars; they also require extremely (almost unrealistically) efficient acceleration 
of relativistic protons to the highest energies within the inner parts of the outflow, and the 
jet kinetic power exceeding the Eddington luminosity by orders of magnitude. 

• The SSC model, which can account for the hard X-ray spectra and also for the synchrotron 
and IC peak frequencies observed in luminous blazars, requires significant departures from 
the minimum power condition in order to explain large luminosity ratio of the low- and 
high-energy spectral components; furthermore, in a dense radiative environment of quasar 
nuclei, energy density of the external radiation fields (provided inevitably by BLR and HDR) 
strongly dominates (in a jet rest frame) over the energy density of the internal jet synchrotron 
photons, and therefore the external- Compton emission is expected to dominate over the SSC 
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emission. 

• The EC models - both EC (BLR) and EC(HDR) variants - can easily account for the large- 
q states of luminous blazars, their hard X-ray spectra, and the observed ratio of synchrotron 
and inverse-Compton peak frequencies; in addition, the main jet parameters implied by 
these models (such as F, Lj, etc.) are consistent with a number of different observational 
constraints. 

• Lack of the bulk-Compton and the Klein-Nishina features in the broad-band spectra of 
luminous blazars, as well as the absence of the low-energy cut-off in their X-ray continua, 
seem to favor the EC(HDR) model over the EC(BLR) variant; however, these particular 
constraints have to be taken with caution, and are not definitive as yet. 

• The most promising observational discrimination between the EC(BLR) and EC(HDR) 
models may be provided by confirmation of the characteristic variability timescale of lumi- 
nous blazars being days/weeks, and by the detection of high-amplitude variability in mm- 
band systematically accompanying 7-ray flares. The previous observational results seem to 
be in fact in agreement with the predictions of the EC(HDR) model; more high quality, well- 
sampled data are however needed in this context, and these are expected to be provided in 
a near future by recently launched Fermi Gamma-ray Space Telescope surveying the whole 
sky, in conjunction with well-sampled multi-band observing campaigns. 

• The electron energy distribution 'reconstructed' from the broad-band emission spectra of 
luminous blazars indicates that the appropriate acceleration processes arc related to rcla- 
tivistic shocks formed in a proton-dominated outflow. This assures self- consistency of the 
favored EC models, in which the blazar emission zone is located far away from the nucleus 
(r > 10^it!g), where the jet may be considered as already fully formed (i.e., accelerated and 
coUimated), and converted to a matter dominated structure. 
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Table 1. Luminous Blazar Sources with the Hardest Recorded X-ray Spectra 



name 
(1) 


z 

(2) 


(3) 


(4) 


(5) 


Ref 
(6) 


S5 0212+73 


2.367 


0.32 + 0.19 






Sambruna et al. (2007) 


PKS 0229+13 


2.059 


0.39 ± 0.09 






Marshall et al. (2005) 


PKS 0413-21 


0.808 


0.39 + 0.12 






Marshall et al. (2005) 


PKS 0528+134 


2.060 


0.12 + 0.26 


1.46 ± 0.04 


1.54 + 0.09 


Donate et al. (2005) 


PKS 0537-286 


3.104 


0.27 + 0.02 


1.47 ± 0.60 




Reeves et al. (2001) 


PKS 0745+241 


0.409 


0.35 + 0.12 






Marshall et al. (2005) 


SWIFT J0746. 3+2548 


2.979 


0.17 + 0.01 






Watanabe et al. (2009) 


PKS 0805-07 


1.837 


0.20 + 0.20 


1.34 + 0.29(7) 




Giommi et al. (2007) 


S5 0836+710 


2.172 


0.34 + 0.04 


1.62 + 0.16 




Donato et al. (2005) 


RGB J0909+039 


3.200 


0.26 + 0.12 






Giommi et al. (2002) 


PKS 1127-145 


1.184 


0.20 + 0.03 


1.70 + 0.31 


1.69 + 0.18 


Siemiginowska et al. (2008) 


PKS 1424-41 


1.522 


0.20 + 0.30 


1.13 + 0.21 




Giommi et al. (2007) 


GB 1428+4217 


4.715 


0.29 + 0.05 






Fabian et al. (1998) 


PKS 1510-089 


0.360 


0.23 + 0.01 


1.47 + 0.21 


1.48 + 0.05 


Kataoka et al. (2008) 


PKS 1830-211 


2.507 


0.09 + 0.05 


1.59 + 0.13 




De Rosa et al. (2005) 


PKS 2149-306 


2.345 


0.38 + 0.08 






Donato et al. (2005) 


PKS 2223+210 


1.959 


0.31 + 0.26 






Donato ct al. (2005) 


3C 454.3 


0.859 


0.34 + 0.06 


1.21 ± 0.06 


1.41 + 0.02 


Donato et al. (2005) 



Note. — (1) Name of a source; (2) Redshift of a source z\ (3) X-ray spectral index a.x\ (4) EGRET 7- 
ray spectral index (Hartman et al. 1999); (5) FERMI 7-ray spectral index (Abdo et al. 2009b); (6) 
References. 
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Fig. 1. — Broadband spectra of fiducial blazars showing the Klein-Nishina effects for the 
EC(BLR) model. Sohd lines are calculated for a power-law type external radiation with 
photon index a — 0, while dashed lines show models with monochromatic external radiation 
(approximated by black-body type emission). Two families of models are shown, correspond- 
ing to the injection of the electron energy distribution oc 7"^-^ with maximal electron 
energy •jmax = 10^, and 10^ (lower and upper curves, respectively). All models are calculated 
for q — 10. Shaded area indicates the energy range observable by the Fermi/LAT. 
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Fig. 2. — Same as Figured! but with q = 100. 
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Fig. 3. — Energy densities of broad line emission {solid line), hot dust radiation {dashed 
line), magnetic field {dotted line), synchrotron radiation {dot-dashed line) and accretion 
disk emission {thin double-dot-dashed line) as seen in the jet comoving frame, as a function 
of distance from the central engine. The shortest timescale of flares corresponding to the jet 
radius i? = r/F at a given distance from the center r, as measured in the observer frame, 
is shown on the top axis. The characteristic radii of broad line region Tblr and dusty torus 
thdr are indicated with a solid and dashed vertical line, respectively. We assumed here: 
jet bulk Lorentz factor F = 20, accretion disk luminosity Ldisk — 10^^ erg/s, synchrotron 
luminosity Lgyn — 10^^ erg/s, magnetic flux Lb — 10^^ erg/s, broad fine region covering 
factor ^BLR — 0.1, dusty torus covering factor ^hdr = 0.1 and dust temperature Td^st — 
10^ K. We assume that beyond vblr stratification of the broad line emission takes form 
dLsLR/dlnr cx 1/r, and hence u'blji oc for r > tblr- Radiation energy density from the 
accretion disk is calculated using formula -u'^ ~ 0.28 {Ldisk/ '^'^^'^c) {Rg/f')'^'^ for Rg of the fast 
rotating black hole with mass Mbh = lO^M© (Dermer & Schlickeiser 2002). 



